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Abstract 

The coupled nonlinear coordinate and spin dynamics of an electron in a double quantum dot with 
spin-orbit interaction is studied semiclassically. The system is driven by an electric field with the 
frequency matching the orbital or the Zeeman resonance in magnetic field. Calculated evolution of 
the spin state is crucially sensitive to the irregularities in the spatial motion and the geometry of 
the host nanostructure. The resulting spin-fiip Rabi frequency has an unusual dependence on the 
field amplitude, demonstrating approach to the chaotic spin motion. In turn, the orbital dynamics 
depends strongly on the spin evolution due to the spin-dependent term in the electron velocity. 

PACS numbers: 72.25.Dc,72.25.Pn,73.63.Kv 
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I. INTRODUCTION 



Spins of electrons confined in nanostructures show a rich coherent and dissipative dy- 
namics (for example,-i^i^) which should be taken into account for the fundamental under- 
standing of their properties and possible application in quantum information^ and spin- 
dependent charge transport. The spin-orbit (SO) coupling, which usually manifests itself 
in low-dimensional systems as the linear in the electron momentum effective field acting 
on electron spin, is one of the main causes of this richness. This momentum-dependent 
coupling offers a new way of manipulating the electron spin by influencing the electron mo- 
mentum, that is by an external periodic electric field. This effect, the electric dipole spin 
resonance,- is in the basis of the very promising proposal of Rashba and Efros^ for spin 
manipulation of electrons confined by and external potential in lateral quantum dots (QDs) 
on the spatial scale of the order between 10 and 100 nm. The weak coupling of the electric 
field to the spin is sufficient to manipulate spins of confined in QDs electrons and holes- and 
produce inhomogeneous spin patterns in lateral semiconductor structures.- The efficiency of 
this technique was demonstrated experimentally with the single-electron QDs.-i^ Another 
side of this coupling is the influence of the spin state and dynamics on the charge motion. 
This effect, studied experimentally and theoretically (see, for example^^^^^) is important, 
especially when it is accumulated at long times leading to qualitatively new transport phe- 
nomena or transitions of the system to well-separated final states. The feedback SO coupling 
effect in the electron weak localization was considered Relaxation of spin states due 

to the coupling to the environmentp^ii^. strongly limits the abilities to initialize, keep, and 
manipulate the desired spin states. 

A fast spin manipulation, often on the time scales strongly limited by the spin relaxation 
times, requires a relatively large electric field amplitude. This electric field, in turn, can 
strongly influence the orbital dynamics and drive the system far beyond the linear motion 
in the confining potential. 

A single QD can be well described by a two-dimensional (2D) parabolic or elliptic poten- 
tial, however, considerably more complex systems are required for applications, such as, at 
least, double quantum dots with a double-minimum potentials. Taking into account the in- 
terest in spin systems in electric field, the full coupled dynamics of the spin and coordinate 
degrees of freedom in a driven QD deserves analysis and understanding of the operation 
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regimes. One of the central questions arising in such a dynamical problem for multi-minima 
potentials is the irregularity and possible transition to the classical and quantum chaos not 
only for the coordinatei^ii^>22 but also for spin. The problems which has to be addressed 
are: how possible irregularities and chaos in the orbital motion influence the spin dynamics 
and what effects in the charge dynamics can be caused by the corresponding irregular spin 
motion ? The analysis of these effects aims on understanding the limits of the spin manipu- 
lation in nanostructures and set a link between spintronics and the nonlinear dynamics and 
chaos theories. 

Most of the studies of spin dynamics in systems not driven by external field have been 
performed so far for the two-dimensional lateral QDs. However, one-dimensional systems, 
including quantum wires and wire-based QDs are of interest^i^i^, promising for the ap- 
plications and can be treated by a thorough analysis of various double-QD problems.— In 
this paper we concentrate on these systems, which demonstrate all important for the under- 
standing features of the coupled spin-charge dynamics in driven systems. 



II. SEMICLASSICAL MODEL FOR COUPLED DYNAMICS 

As a model we consider a particle with mass m in a quartic potential 

u^,)^u„[-2{^\{^y (1) 

describing a double-minimum system where the minima are separated by 2d in space with 
the barrier of Uq in energy. Below we refer to this system as to the double QD (DQD) 
system with each region U{x) < considered as a one-dimensional QD. For this potential 
the classical^ and quantum tunneling-determined^^ dynamics driven by electric field have 
been studied. Here we concentrate on the classical dynamics by considering mainly a wide 
structure with d = 100-\/2 nm and Uq = 25 meV. This high energy scale is on the order 
of 300 K, and, therefore, at temperatures below 100 K the orbital dynamics is only weakly 
sensitive to the temperature. For this reason, the thermal effects on the electron momentum 
and kinetic energy, including the spread and the activated over-the-barrier motion, will be 
neglected. The frequency of the single QD small-amplitude oscillations Uq = 2\^^Uo/m/d 
is 4.96- 10^^ s~^ at GaAs electron effective mass m = 0.067mo (mo is the free electron mass), 
the semiclassically calculated number of energy levels in a single QD Ne<o = 4v^(i-\/mt/o/37r 
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is close to eight, and the interdot semiclassical tunnehng probabihty exp 

is vanishingly small, which allows to treat the dynamics semiclassically (we put h = 1). 

To characterize the SO interaction, we use the Dresselhaus type of coupling present in 
all the zincblend structures in the form: 

^so = OiD (i^'O-) , (2) 

where ao is the bulk Dresselhaus constant, cr are the Pauli matrices, = Px{v1 ~ Pz) 
with other components obtained by the cyclic permutation, and pi are the corresponding 
momentum components: = —id/dx, etc. In a one-dimensional structure, extended along 
the X— axis, where the transverse motion is quantized, the Dresselhaus term is reduced to 
the linear in momentum term jSa^Px, where (3 = aniPy ~ pD with the quantum mechanical 
expectation value {p^ — pi) for the ground state of the quantized transverse motion. — The 
value of (3 depends on the asymmetry of the system cross-section and vanishes if it has a 
perfect square or circle shape. In the GaAs-based structures with the transverse dimensions 
of few nanometers, (3 is on the order of 0.1-0.5x10^^ eVcm. 

In the presence of the 2;-axis oriented magnetic field Bz the ID Hamiltonian acquires the 
Zeeman term only, and with the electric field E{t) applied along the OC clXlS it has the form 

H=^ + U{x) - eE{t)x + (3a^px + ^^b^^^.. (3) 
Im 2 

Following the approach of Refj6|, the direct Zeeman coupling of the magnetic field of the 
electromagnetic wave to the electron spin was neglected in Eq.([3]). The evolution of mean 
values for variable X is governed by the equation of motion X = i[H,X] where [, ] stands 
for the commutator. By applying the commutation rules for x, p^ and for the Pauli matrices 
one obtains the following coupled equations for the coordinate and spin dynamics: 

X = Px/m + /?cr^, 

Px = -dU /dx + eE{t) - p^/Tp, 

& = [(t^soX + sgn{g)ujLz) x cr] , (4) 

where iVso = 2/3px, uj^ = IglfisB^ is the Larmor frequency, and x, z are the unit vectors, 
and upper dot stands for the time derivative. Since in the chosen range of parameters, the 
orbital {x,px) dynamic can be treated classically, here the momentum and coordinate are 
considered as the quantities governed by the classical equations of motion. We include here 



the momentum relaxation with Tp = 5 ps being a typical relaxation time for the conduction 
electron. Since oooTp ^ 1, this relaxation influences the dynamics weakly. Also, we neglect 
the spin relaxation which occurs on the time scales^^ much longer than the maximum times 
considered in this paper. Hence, the spin states studied here can be considered as long- 
standing spin states which properties are of particular importance for spintronics. The 
system (j4]) should be accompanied by the initial conditions for the position x(0), momentum 
Px{0), and spin components cr*(0) which values are determined in real structures by the 
method of the spin state preparation before the driving field is applied. The typical set of 
the initial parameters which we consider in the paper describes the spin-up electron located 
at the bottom of a particular quantum well. 

III. DYNAMICS IN ELECTRIC FIELD 

A. Zeeman resonance: uj = ul 

To begin with we consider analytically the resonance with E{t) = —Eq sin{uJi,t) at small 
Eq. For the moment, we suppose that \e\EQ ^ Uo/d, such that all nonlinear terms in the 
force acting on the electron are small compared to the linear ones and the orbital oscillations 
x{t) are very close to the harmonic. Near the spin-flip resonance, we can consider the spin 
dynamics with the explicitly spin- dependent part of the Hamiltonian: 



in the two-level model, taking only the resonant contribution exp{—iujLt) in sm{ujLt) = 
{exp{iujLt) — exp(— ztu^t)) /2i. In this Rotating Wave Approximation neglecting contribu- 
tions of the off-resonance terms, the Hamiltonian ([5]) can be written as: 



with the Rabi frequency Q = Ppr where pr = IcEqIul/uq is the amplitude of the driven 
momentum. The solution to the Schroedinger equation with Hamiltonian ([6]) yields 

(cr^) = — sin(nt) sinuit, (a^) = — sin(fit) cosc^it. 



(5) 



s 






(7) 
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The spin expectation values ([7j) show the Rabi oscillations in while the amplitude of the 
and Larmor precession at ujl is modulated with the Rabi frequency. 
This analytical result can be compared to the numerical solution of system (jlj) with initial 
conditions corresponding to the spin-up electron in the left QD. Examples of dynamics 
and Fourier power spectra are shown in Fig{T] for g = —0.45 (we assume the bulk GaAs 
value here), f3 = 0.2 ■ 10~^ eVcm, and 5^ = 4 T corresponding to ujl = 1.6 ■ 10^^ s~^. The 
power spectra v^j^ were obtained as the Fourier transform of the numerical solution for the 
dynamical variables^^ u(t = n = 1, . . . , of (jl]) at the frequencies Um = 27rm/tAr defined 
by the overall time ^at of the observation: 

N 

v.^ = Y.<in)e-'-'^--'^^'^~'^"'. (8) 

ra=l 

The time dependence is shown in Fig{T] (a) for the exact Zeeman resonance uj = ujl for 
three cases illustrating different regimes of the dynamics, the role of the structure geometry, 
and the effect of the interdot transfer for the spin dynamics and the Rabi frequency. 

When the electron is still confined to one dot, the spin motion is regular and the Rabi 
frequency is approximately linear in the external field amplitude. When the amplitude is 
slightly increased from 2.55 kV/cm to 2.6 kV/cm, the particle can travel through both QDs 
causing a strong change in the Rabi frequency and the entire spin dynamics. The Rabi 
frequency in this case is not related directly to ujq or the geometrically defined SO frequency 
ujg = jSpg {pg = \/ 2mU q) sluce the evolution here is strongly nonlinear. Above the transfer 
threshold the motion is the superposition of two processes shown in Fig. (2): a slow interdot 
oscillations with the corresponding frequency uol and fast oscillations at the frequencies on 
the order of ujq ^ ujl in the vicinity of the minima, not leading to the spin resonance. The 
increase in the maximum of the momentum compared the lower-field case, where the 
transfer is still prohibited and the motion occurs at the frequencies determined by ujl is on 
the order of uoq/ujl. However, since this very large increase corresponds to the frequencies 
far away from the Zeeman resonance, it influences the spin dynamics relatively weakly. The 
momentum Fourier component corresponding to the slow motion is proportional to 
now, causing the abrupt change in the Rabi frequency. The spin behavior becomes strongly 
irregular. In FigHl^b) we present the corresponding Fourier power spectra in the case of in- 
terminimum transfer showing this spin dynamics. If a broader DQD structure is considered, 
then a lower driving field can excite the oscillations with larger electron momentum due to 
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FIG. 1: (Color online) (a) Evolution of cr^(t) at uj = u;l for confined motion in the left QD at 
£"0 = 2.55 kV/cm (solid black line) and £"0 = 2.6 kV/cm when the transfer is induced (line marked 
n, red) for the d = 100^2 nm structure, and for a a structure with d = 200\/2 nm with the transfer 
at smaller Eq = 1.5 kV/cm (line marked w, blue). Dashed black line shows cr^{t) at Eq = 2.55 
kV/cm, d = WOV2 nm, UJ = l.lujL. The full spin flip is achieved only at the exact resonance, (b) 
Fourier transform spectra for Eq = 2.6 kV/cm, d = 100\/2 nm (line marked n, red) and Eq = 1.5 
kV/cm, d = 200\/2 nm (line marked w, blue) showing the approach to the irregular spin dynamics, 
(c) Phase portrait for uj = ujl, d = 100\/2 nm, Eq = 2.55 kV/cm, time t < 700 ps, (d) same for 
Eq = 2.6 kV/cm. 
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FIG. 2: (Color online) Phase portrait {x,px) of the system at lo = lol and Eq = 2.6 kV/cm when 
the transfer is possible at time t < 700 ps. The pattern is dominated by the motion in a single 
(left or right) QD at the frequencies on the order of ujq. The interdot transfer is a relatively rare 
event at the frequency lol- 

a more extended and smoother confining potential. Here the Rabi frequency of cr^ flip is 
increased by a factor of two when comparing the d = 100^/2 and d = 200^2 nm structures 
since the interwell transfer is characterized by a factor of two greater amplitude of px at ul. 
The corresponding phase portraits are presented in Figll](c) and Figll](d) for the regular and 
irregular dynamics shown in Fig.l. The plot in FiglT](d) is a signature of approach to chaotic 
regime, where the back and forth spin motion on a short time scale is superimposed by a 
more regular dynamics leading to the spin-flip on the longer time scale. It should be noted 
that the terms "chaotic motion" together with the "irregular motion" used in the paper 
refer to the visual shapes of the parameter dynamics rather than to a strictly defined dy- 
namical chaos. The possible appearance and investigation of dynamical chaos in our system 
is outside of the scope of the present paper since our goal here is limited to observation and 
description of the transition to irregular dynamics for coupled coordinate and spin degrees of 
freedom. Hence, we do not perform a unnecessarily detailed investigation of the dynamical 
properties including, for example, the analysis of Lyapunov exponent spectrum.— Still, the 
analysis presented above and in FigH] certainly lead to a conclusion that the spin dynamics 
in our system is strongly affected by the SO coupling and can be tuned by varying both the 
geometrical parameters of the nanostructure and the amplitude of external electric field. 
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B. Orbital resonance: tu = uq 



Another type of dynamics to be considered is the evolution under the orbital resonance 
field at the QD oscillation frequency c^o- When the electric field amplitude is high enough 
to overcome the barrier and the momentum relaxation, the particle can be driven through 
both QDs which is shown in Figl3](a) for Eq = 1.6 kV/cm, although a full fiip for is 
not achieved even if an artificially large g = 4.0 is considered (dashed lines) as shown in 
FiglSj^b). The corresponding Fourier power spectra for g = —0.45 are shown in FigjSj^c). 
By comparing Fig|3] and Fig{T] it becomes clear that the expansion in the space motion 
leads to the enrichment in both coordinate and spin power spectra due to SO coupling. The 
dominating part of spectrum is located between the lowest ul and the highest Uq frequencies 
of the system which gives rise to various contributions in the whole interval between them. 
The phase portrait for (cr^(t), cr^(t)) is presented in Figj3]^d) where one can see that the 
spin dynamics indeed does not show a full spin fiip and is asymmetric with respect to the 
maximum spin projections since the potential U{x) is spatially asymmetric with respect to 
the single QD minimum. 

The character of the spin evolution changes significantly in a broader DQD structure. Let 
us consider d = 200^2 nm which is twice the initial value and the corresponding frequency 
u)q is now factor of two lower. By a moderate increase in Eq from 1.6 to 2.0 kV/cm one can 
significantly modify the particle dynamics and, in particular, achieve a full fiip shown 
in Figl3](e). The calculations show that after some period of transient dynamics a steady 

oscillations are reached which is further confirmed by the Fourier power spectrum the 
dominating peak at new ujq. The corresponding phase portrait for (cr^(t), cr^(t)) with a full 
spin fiip is shown in Figl3](e) which comparison to Fig{T](d) indicates that the spin evolution 
here has a more regular character. This can be also confirmed by the analysis of cr^(t) 
dynamics which shows regular oscillations after the transient period. The results presented 
above demonstrate again that the evolution of a spin in the presence of SO interaction is 
very sensitive to the spatial motion and geometrical parameters of the host nanostructure. 

To study the long-term effect of the orbital resonance field, a finite length pulse with 



will be considered. We focus on the asymptotic values at t tp for x(t) and <J^(t) for the 
coupled spin-coordinate evolution. Thus, we take to = 2tp and tp varying from 2To to 12To, 




(9) 
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FIG. 3: (Color online) Evolution of (a) x{t) and (b) cr^(t) at the primary resonance uj = loq and 
Eq = 1.6 kV/cm where the particle with g = —0.45 (solid black line) and g = 4.0 (dashed red line) 
is driven through both QDs. The increase in the (7-factor leads to much stronger modulation, 
(c) The Fourier power spectra corresponding to x{t) (dashed red line) and cr^{t) (solid blue line) 
at g = —0.45 shown in (a) and (b), respectively, (d) Phase portrait for u = loq, d = 100\/2 nm, 
Eq = 1.6 kV/cm, time t < 50 ps (e) Same for d = 200\/2 nm, Eq = 2.0 kV/cm, demonstrating the 
achieved full spin flip. 
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FIG. 4: (Color online) Asymptotic values xqd = x{t oo)/d (filled circles) and a^{t — > oo) (filled 
squares) after the system has been driven by resonant pulses with Eq = 3 kV/cm and varying 
duration for (a) = and (b) = 4T. The magnetic field strongly affects both the spin and 
the coordinate dynamics demonstrating the role of SO coupling. 

where To = 27t/ujq, and the amplitude Eq = 3 kV/cm. The results are shown in FigJD^a) 
for the zero magnetic field and in FigHl^b) for Bz = 4T. It follows from (jlj) that at zero 
magnetic field = and for the given initial condition x(0) = the solutions for x{t) and 
a^(t),(T^(t) are related as 



ay 






a' 


» _ 





cos {2f3mx) — sin {2(3mx) 
sin {2(3mx) cos (2/3mx) 



a^(0) 
a^(0) 



(10) 



Therefore, the spin returns to the initial state along with the coordinate as it can be seen in 
FigSl^a). When the magnetic field is applied, the spin-orbit jSa^Px and Zeeman {g/2)fiB<^^Bz 
terms in Eq.([5]) do not commute, and such relation is no longer valid. Here, both the final 
x{t) and cr^{t) values form a different pattern shown in FigJl](b). We emphasize that not only 
the spin but the coordinate asymptotic also becomes magnetic field-dependent due to the 
SO coupling term f3a^ in the velocity, as it can be seen by comparing FiglD^a) and (b). This 
result is an important complementary effect of the influence of spin dynamics on the charge 
evolution which should be taken into account when both the charge and spin dynamics in 
systems with strong spin-orbit coupling are of interest. 
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IV. CONCLUSIONS 



We have studied in the semiclassical approximation coupled coordinate-spin dynamics of 
electron in a one-dimensional double quantum dot with the Dresselhaus type of SO coupling 
in external magnetic field. The considered system was driven by a harmonic electric field 
with the frequency matching the transition between the orbital levels of the quantum dot 
or the Zeeman resonance, or by an intense finite-time orbital resonant field pulse. The 
increase in the field amplitude makes the electron transfer between the potential minima 
possible and triggers strongly irregular behavior in both, coupled charge and spin, channels. 
These irregularities in the spin motion suggest that in a nonlinear system evolution of spin 
coupled to the momentum is very sensitive to the spatial motion. The spin flips, therefore, 
cannot be well described by a single Rabi frequency. In turn, the anomalous spin- dependent 
contribution to the electron velocity leads to the position time dependence x{t) strongly 
different from that expected for the zero spin-orbit coupling. The study of the development 
of these irregularities into the real chaos in the spin subsystem, involving the nontrivial 
Lyapunov spectrum analysis, is an interesting problem for the future research. The geometry 
of the nanostructure strongly influences the spin dynamics in this nonlinear regime. These 
conclusions emphasize the importance of the possible nonlinear behavior for the initialization 
and manipulation of spin states in quantum dots in the experimentally accessible range of 
system parameters. 
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